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Abstract 

We perform a detailed analysis of lepton flavour violation (LFV) within minimal 
see-saw type extensions of the Standard Model (SM), which give a viable mechanism 
of neutrino mass generation and provide new particle content at the electroweak scale. 
We focus, mainly, on predictions and constraints set on each scenario from [i — > ej, 
[i — > 3e and /i — e conversion in the nuclei. In this class of models, the flavour structure 
of the Yukawa couplings between the additional scalar and fermion representations 
and the SM leptons is highly constrained by neutrino oscillation measurements. In 
particular, we show that in some regions of the parameters space of type I and type 
II see-saw models, the Dirac and Majorana phases of the neutrino mixing matrix, the 
ordering and hierarchy of the active neutrino mass spectrum as well as the value of the 
reactor mixing angle #13 may considerably affect the size of the LFV observables. The 
interplay of the latter clearly allows to discriminate among the different low energy 
see-saw possibilities. 



1 Introduction 

After several decades of neutrino experiments, a clear quantitative picture of the neutrino 
oscillation parameters is gradually emerging (see, e.g. [1]). The Super-Kamiokande col- 
laboration established that the atmospheric neutrino mass squared splitting is |Am^| ~ 
O(10~ 3 eV 2 ) and that the corresponding mixing angle is large, possibly maximal #23 — 7r /4 
[2]. The data from SNO, Super-Kamiokande and KamLAND experiments [31 IH |5] allowed 
to established the large mixing angle solution as a unique solution of the long standing so- 
lar neutrino problem, with a solar neutrino mass squared splitting AttIq ~ O(10~ 5 eV 2 ) 

1 Also at: Institute of Nuclear Research and Nuclear Energy, Bulgarian Academy of Sciences, 1784 Sofia, 
Bulgaria 
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and mixing angle #12 — arcsin(-y/0.3). A series of subsequent experiments, using reactor 
and accelerator neutrinos, have pinned down the atmospheric and solar neutrino oscillation 
parameters with a few to several percent accuracy, as summarized in Table [TJ 

Furthermore, in June of 2011 the T2K collaboration reported [Hj evidence at 2.5a for a 
non-zero value of the angle #13. Subsequently the MINOS [7] and Double Chooz [5] collab- 
orations also reported evidence for # 13 7^ 0, although with a smaller statistical significance. 
Global analyses of the neutrino oscillation data, including the data from the T2K and MI- 
NOS experiments, performed in [9], [TO] showed that actually sin0 13 7^ at > 3a. The results 
of the analysis (9], in which Am^ = AttIq and jAm^j = |Am|J were determined as well, 
are shown in Table [TJ 

Recently, the first data of the Daya Bay reactor antineutrino experiment on 6*13 was 
published [TTJ. The value of sin 2 2# 13 was measured with a rather high precision and was 
found to be different from zero at 5.2a: 

sin 2 29 13 = 0.092 ± 0.016 ± 0.005 , 0.04 < sin 2 29 13 < 0.14 , 3a, (1) 

where we have given also the 3a interval of allowed values of sin 2 29i 3 . Subsequently, the 
RENO experiment reported a 4.9<r evidence for a non-zero value of #13 |12j : 

sin 2 29 13 = 0.113 ±0.013 ±0.019. (2) 

The value of #13 determined in the RENO experiment is compatible with that measured in 
the Daya Bay experiment. It is interesting to note also that the mean value of sin 2 # 13 found 
in the global analysis of the neutrino oscillation data in [9] differs very little from the mean 
values found in the Daya Bay and RENO experiments. 

The results on #13 described above will have far reaching implications for the program of 
research in neutrino physics. A relatively large value of sin 6*13 ~ 0.15 opens up the possibili- 
ties, in particular, i) for searching for CP violation effects in neutrino oscillations experiments 
with high intensity accelerator neutrino beams (like T2K, NOi^A, etc.); ii) for determining 
the sign of Am 2 2 , and thus the type of neutrino mass spectrum, which can be with normal or 
inverted ordering (see, e.g. [I]), in the long baseline neutrino oscillation experiments at accel- 
erators (NOi/A, etc.), in the experiments studying the oscillations of atmospheric neutrinos 
(see, e.g. [U]), as well as in experiments with reactor antineutrinos [JBj. It has important 
implications for the neutrinoless double beta ((/3/3) 0i y-) decay phenomenology in the case of 
neutrino mass spectrum with normal ordering (NO) [IT]. A value of sin #13 ^ 0.09 is a nec- 
essary condition for a successful "flavoured" leptogenesis when the CP violation required for 
the generation of the matter-antimatter asymmetry of the Universe is provided entirely by 
the Dirac CP violating phase in the neutrino mixing matrix [15]. I 2 ! As was already discussed 
to some extent in the literature and we will see further, in certain specific cases a value 
of sin 0x3 ~ 0.15 can have important implications also for the phenomenology of the lepton 
flavour violation (LFV) processes involving the charged leptons in theories incorporating one 
of the possible TeV scale see-saw mechanisms of neutrino mass generation. 

Despite the compelling evidence for the nonconservation of the leptonic flavour in neutrino 
oscillations, all searches for lepton flavour violation (LFV) in the charged lepton sector have 

2 If indeed sin 13 = 0.15 and the neutrino mass spectrum is with inverted ordering (IO), further important 
implications for "flavoured" leptogenesis are possible |19j . 
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Table 1: The best-fit values and 3a allowed ranges of the 3-neutrino oscillation parameters, 
derived from a global fit of the current neutrino oscillation data, including the T2K and 
MINOS (but not the Daya Bay) results (from [5]). The Daya Bay data [TT] on sin 2 # 13 is 
given in the last line. The values (values in brackets) of sin 2 # 12 are obtained using the "old" 
[T3] ("new" jH]) reactor v e fluxes in the analysis. 

Parameter best-fit (±lo~) 3cr 

Am 2 Q t 10 " 5 ey2 ] 7.58l°; 22 6.99 - 8.18 

\Am 2 A \ [KT 3 eV 2 ] 2.35^ 2.06 - 2.67 

sin 2 # 12 0.3061°;°^ 0.259(0.265) - 0.359(0.364) 

sin 2 9 23 0.42±g;g| 0.34-0.64 

sin 2 # 13 [9] 0.021(0.025) ±g:gg? 0.001(0.005) - 0.044(0.050) 

sin 2 0i3 III] 0.0236 ±0.0042 0.010 - 0.036 



been so far negative. The best limits follow from the non-observation of the LFV muon 

decays /i + — » e + 7 and /i + — > e + e~e + , 

BR(/i + -»■ e + 7 ) < 2.4 x 10~ 12 [20] , (3) 
BR(/i+ -»■ e + e"e + ) < 1.0 x 10~ 12 [21] , (4) 

and from the non-observation of conversion of muons into electrons in Titanium, 

CR(/xTi -»■ eTi) < 4.3 x 10~ 12 [22] . (5) 

Besides, there are stringent constraints on the tau-muon and tau-electron flavour violation 
from the non-observation of LFV radiative tau decays |23j : 

BR(r — »■ fWj) < 4.4 x 10" 8 , (6) 
BR(r ->■ erf) < 3.3 x 10" 8 . (7) 

The indicated stringent upper limits on the rates of the LFV processes involving the 
charged leptons lead to severe constraints on models of new physics which predict new 
particles at the electroweak scale coupled to the charged leptons. Indeed, the dipole operator 
which leads to the process /i — > e^y has the form: 

- £ = mM^ + l,t E \)a vp eF yp + h.c. (8) 

where and are, respectively, the transition magnetic and electric dipole moment 
form factors. This operator is generated at the quantum level through particles with masses 
A which couple to the charged leptons, hence the form factors can be parameterised as 
= 16 t A 2 , Ofie being a parameter which measures the strength of the coupling of the new 
particles to the electron and the muon. The present experimental limit on BR(/i — > ej) sets 
the following upper limit on the two form factors: |/^|, \ Jmi\ ~ 10~ 12 GeV -2 . The latter in 
turn translates into A > 20 TeV if 6^ ~ 1/^2, or in 6^ < 0.01 if A ~ 300 GeV. It is then 
apparent that experiments searching for lepton flavour violation can probe models of new 
physics which cannot be tested in collider experiments, either because the new particles are 
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not kinematically accessible with the available collider energies, or because the couplings of 
the new particles to the Standard Model (SM) particles are too feeble to produce the former 
with rates necessary for their observation given the luminosity of the present colliders. 

Low scale see-saw models are a particularly interesting class of models of new physics 
which are severely constrained by experiments searching for lepton flavour violation. In this 
class of models the flavour structure of the couplings of the new particles to the charged 
leptons is basically determined by the requirement of reproducing the data on the neutrino 
oscillation parameters [2H (23 [261 123 [28]. Hence, the rates for the lepton flavour violating 
processes in the charged lepton sector can be calculated in terms of a few parameters, the 
predicted rates being possibly within the reach of the future experiments searching for lepton 
flavour violation, even when the parameters of the model do not allow production of the new 
particles with observable rates at the LHC [28] . 

The role of the experiments searching for lepton flavour violation to constrain low scale 
see-saw models will be significantly strengthened in the next years. Searches for \i — e 
conversion at the planned COMET experiment at KEK [29] and Mu2e experiment at Fer- 
milab [3D] aim to reach sensitivity to CR(/iAl — > eAl) rs 10~ 16 , while, in the longer run, the 
PRISM/PRIME experiment in KEK [31] and the project-X experiment in Fermilab [32] are 
being designed to probe values of the \i — e conversion rate on Ti, which are by 2 orders of 
magnitude smaller, CR(/iTi — > eTi) w 10 -18 [31]. If these experiments reach the projected 
sensitivity without observing a signal, the upper limits on the form factors /mi, Jei will 
improve by two orders of magnitude. There are also plans to perform a new search for the 
/i + — y e + e~e + decay [33], which will probe values of the corresponding branching ratio down 
to BR(/i + — > e + e~e + ) ~ 10~ 15 , i.e., by 3 orders of magnitude smaller than the best current 
upper limit eq. (jl]). Furthermore, searches for tau lepton flavour violation at superB factories 
aim to reach a sensitivity to BR(r — > (/i, e)j) ~ 1CT 9 [3U 155] . 

In this paper we will study the constraints on low (TeV) scale see-saw models of neutrino 
mass generation from present and future experiments searching for lepton flavour violation, 
with especial emphasis on /i — e conversion in nuclei, which is among all search strategies 
the one with brightest perspectives. The paper is organized as follows: in Sections El [3] and 
H] we review the main features of the three types of see-saw mechanisms. We discuss for 
each scenario the predictions and experimental constraints on the relevant parameter space 
arising from LFV processes. The results are summarized and discussed in the concluding 
Section [5j 



2 TeV Scale Type I See-Saw Model 

We consider in detail in this Section LFV processes emerging in type I see-saw extensions of 
the SM [36] . We denote the light and heavy Majorana mass eigenstates of the type I see-saw 
model as Xi an d Nk, respectively. El The charged and neutral current weak interactions 
involving the light Majorana neutrinos have the form: 

C-CC = ~ A= ha VlL W a + h.C. = " 4| ha ((1 + vW) ti XiL W" + h.C. , (9) 

Cnc = -^-T^laVi L Z a =-J-xIEla (W(l + rj + ^)U) XjL Z a , (10) 
3 We use the same notations as in [281 [57] . 
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where (1 +rj)U = Upuns is the Pontecorvo, Maki, Nakagawa, Sakata (PMNS) neutrino mix- 
ing matrix [38| [39] . U is a 3 x 3 unitary matrix which diagonalises the Majorana mass matrix 
of the left-handed (LH) flavour neutrinos vn (generated by the see-saw mechanism), and 
the matrix rj characterises the deviations from unitarity of the PMNS matrix. The elements 
of f/pMNS are determined in experiments studying the oscillations of the flavour neutrinos 
and antineutrinos, vg_ and Pg, £ = e, fi, r, at relatively low energies. In these experiments the 
initial states of the flavour neutrinos, produced in weak processes, are coherent superposi- 
tions of the states of the light massive Majorana neutrino Xi only. The states of the heavy 
Majorana neutrino Nj are not present in the superpositions representing the initial flavour 
neutrino states and this leads to deviations from unitarity of the PMNS matrix. 

The matrix rj can be expressed in terms of a matrix RV whose elements (RV)e k determine 
the strength of the charged current (CC) and neutral current (NC) weak interaction couplings 
of the heavy Majorana neutrinos Nk to the W^-boson and the charged lepton £, and to the 
Z° boson and the left-handed (LH) flavour neutrino z/^, I = e,fi,r: 

C-cc = -^h a (RV)tk(l-7s)N k W a + h.c, (11) 
C N NC = -j-VEla{RV), k {l-i 5 )N k Z a + h.c. (12) 

Here V is the unitary matrix which diagonalises the Majorana mass matrix of the heavy RH 
neutrinos and the matrix R is determined by (see [37]) R* = M^ 1 , Mp and being the 
neutrino Dirac and the RH neutrino Majorana mass matrices, respectively, \Mjj\ <C |Mjy|. 
We have: 

rj = - l -RR) = ~(RV)(RVy = 77+ . (13) 

The elements of the matrices RV and rj can be constrained by using the existing neutrino 
oscillation data, data on electroweak (EW) processes, etc. [HH EH]. They should satisfy also 
the constraint which is characteristic of the type I see-saw mechanism under discussion: 

| J2( RV Ti'k M * ( R V)l\ = I K)«| < 1 eV , £',£ — e,fj,,r . (14) 
k 

Here m v is the Majorana mass matrix of the LH flavour neutrinos generated by the see-saw 
mechanism. The upper limit \(m u )£>e\ < 1 eV, £,£' = e,fi,T, follows from the existing data 
on the neutrino masses and on the neutrino mixing [12]. For the values of the masses M k 
of the heavy Majorana neutrinos N k of interest for the present study, M k = (9(100 — 1000) 
GeV, the simplest scheme in which the constraint (fT4|) can be satisfied is [28J the scheme 
with two heavy Majorana neutrinos (see, e.g., [131 HH 05] ) , N± and N 2 , which form a pseudo- 
Dirac neutrino pair [l6[ H7J: M 2 = Mi(l + z), where 2<1, which naturally arises in type 
I see-saw models with a mildly broken lepton number symmetry [25J and in the inverse see- 
saw model [1SI HH]- In the scenario where the CC and NC couplings of iVx 2 are sizable, the 
requirement of reproducing the correct low energy neutrino oscillation parameters constrains 
significantly [25[ [26] and in certain cases determines the Yukawa couplings [2U [27J [28] . 
Correspondingly, the elements (RV)a and (RV)^ in eqs. (fTTj) and (JT2J) are also determined 
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and in the case of interest take the form [281: 



\(RV) n \ 2 = I V lf 2 ^ U e3 + iy/m 2 /m 3 U e2 ' , NH , (15) 
2 M( m 2 + m 3 



\(RV) a f = ^^^^|^ + ^VW^^r^ |^|C/, 2 + ^x! 2 ,IH, (16) 
(RV) e2 = ±i(RV) n \h^, t = e,n,r, ( 17 ) 



M- 
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where y represents the maximum eigenvalue of the neutrino Yukawa matrix and v ~ 174 
GeV. In eq. ([16]) we have used the fact that for the IH spectrum one has m x = m 2 . Due to 
the relation (TlTl) between (RV)a and (RV)e 2 , eq. (fT4l) is automatically satisfied. 

Upper bounds on the couplings of RH neutrinos with SM particles can be obtained 
from the low energy electroweak precision data on lepton number conserving processes like 
Ti — > £n#, Z y isls and other tree- level processes involving light neutrinos in the final state 
[4"0] . These bounds read: 

\(RV) el \ 2 < 2xl(T 3 , (18) 
\(RV)»i\ 2 < 0.8 xlO" 3 , (19) 
\(RV) Tl \ 2 < 2.6 x 10" 3 . (20) 

Let us add that in the class of type I see-saw models with two heavy Majorana neutri- 
nos we are considering (see, e.g., [I3j HH US]), one of the three light (Majorana) neutrinos 
is massless and hence the neutrino mass spectrum is hierarchical. Two possible types of 
hierarchical spectrum are allowed by the current neutrino data (see, e.g., p]): i) normal 
hierarchical (NH), m\ = 0, m 2 = ^ArriQ and m 3 = ^ Am\ , where Amg = m 2 — m\ > 



and Am\ = m\ — mf; ii) inverted hierarchical (IH), m 3 = 0, m 2 = y/\Am\ \ and 
mi = y?\Am\ | - AttIq = \J\Am\ |, where AttIq = m\—m\ > and Am\ = m\—m\ < 0. 
In both cases we have: Am 2 & /\Am\ \ = 0.03 <C 1. 

The numerical results we will present further will be obtained employing the standard 
parametrisation for the unitary matrix U : 

U = V(9 12 ,9 23 ,9 13 , 5) Q(a 21 ,a 31 ). (21) 

Here (see, e.g., [I]) 

s 13 e~ lS \ ( cia sia 
V = | c 23 s 23 | | 1 II -s 12 c 12 | , (22) 

1 

where we have used the standard notation c^- = cos 9ij, Sy = sin%, 5 is the Dirac CP 
violation (CPV) phase and the matrix Q contains the two Majorana CPV phases Q [50] , 

Q = diag(l,e iQ2l/2 ,e ia3l/2 ). (23) 





In the case of the type II see-saw mechanism, to be discussed in Section 3, we have r\ = and thus the 
neutrino mixing matrix coincides with U: ?7pmns = U. We will employ the parametrisation (f2Tj) - (|23[) for 
U also in that case. 
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We recall that C/pmns = (1 + ifjU. Thus, up to corrections which depend on the elements 
of the matrix 77 whose absolute values, however, do not exceed approximately 5 x 10 -3 
|40j . the values of the angles 612, #23 and #13 coincide with the values of the solar neutrino, 
atmospheric neutrino and the 1-3 (or "reactor" ) mixing angles, determined in the 3- neutrino 
mixing analyses of the neutrino oscillation data and reported in Table [TJ 

Given the neutrino masses and mixing angles, the TeV scale type I see-saw scenario we 
are considering is characterised by four parameters [28]: the mass (scale) M 1; the Yukawa 
coupling y, the parameter z of the splitting between the masses of the two heavy Majorana 
neutrinos and a phase u. The mass Mi and the Yukawa coupling y can be determined, 
in principle, from the measured rates of two lepton flavour violating (LFV) processes, the 
/i — > e7 decay and the fi — e conversion in nuclei, for instance. The mass splitting parameter 
z and the phase u, together with Mi and y, enter, e.g., into the expression for the rate of 
(/3/9)oy-decay, predicted by the model. The latter was discussed in detail in [28] . 



2.1 The 11 — y cy Decay 

In this subsection we update briefly the discussion of the limits on the parameters of the 
TeV scale type I see-saw model, derived in [28] using the experimental upper bound on the 
fi — > ey decay rate obtained in 1999 in the MEGA experiment [5T] . After the publication 
of [28], the MEG collaboration reported a new more stringent upper bound on the fi — > &y 
decay rate [2U] given in eq. ([3]). Such an update is also necessary in view of the relatively large 
nonzero value of the reactor angle #13 measured in the Daya Bay and RENO experiments 
[TT| [T2] and reported in eqs. ([1]) and ([2]). As was discussed in [28], in particular, the rate of 
the fi — > e'-f decay in the type I see-saw scheme considered can be strongly suppressed for 
certain values of 6> 13 . 

The fi — > cy decay branching ratio in the scenario under discussion is given by [52J 153"] : 

BR(^e 7) = r . r("-^> = ^=|71», (24) 
r(/i —^e + u^ + Ue) 327r 

where a cm is the fine structure constant and [2B~| 

2 + z 



1 + 2 



(RV); i (RV) e i\\G(X)-G(0)\ . (25) 



In eqs. ff2^|) and (T25]) the loop integration function G(x) has the form: 

N 10 - 43x + 78a; 2 - 49x 3 + 4x 4 + I82 3 log(x) 
G{X) = 3£=i)i ' (26) 

where X = {Mi/My/) 2 . In deriving the expression for the matrix element T, eq. (125 p . we 
have assumed that the difference between Mi and M 2 is negligibly small and used Mi = M 2 . 
It is easy to verify that G(x) is a monotonic function which takes values in the interval 
[4/3, 10/3], with G(x) ^ f - x for x < 1. 

Using the expressions of | (RV)^i \ 2 and | (RV) e i | 2 in terms of neutrino parameters, eqs. ( |T5l) 
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and (fT6|) . we obtain the /i — > e'y decay branching ratio for the NH and IH spectra: 



NH: BR(/i^e 7 ) 
3a em / y 2 v 2 m 3 



32tt V Mf m 2 + m 3 




Ue3 + i\l—U e2 
m 3 



2 



[G(X) - G(0)} 2 , (27) 



IH : BR(yU e 7 ) 

W (Sf ^) V,2 + ^ M i| 2 |^ e2 + ^ e i| 2 [G(X) - G(0)] 2 . (28) 



The data on the process \i — > cy set very stringent constraints on the TeV scale type I 
see-saw mechanism. The current upper bound on BR(/i — > cy) was obtained in the MEG 
experiment at PSI [20J and is given in eq. ([3]). It is an improvement by a factor of 5 of the 
previous best upper limit of the MEGA experiment, published in 1999 [5.1]. The projected 
sensitivity of the MEG experiment is BR(/i ej) ~ 10" 13 [20]. For M x = 100 GeV (Mi = 1 
TeV) and z < lwe get the following upper limit on the product \ (RV)* (jLl (RV) e \ \ of the heavy 
Majorana neutrino couplings to the muon (electron) and the boson and to the Z° boson 
from the the upper limit eq. ([3]): 

| (RV)*^ (RV) el I < 0.8 x 10" 4 (0.3 x 10~ 4 ) , (29) 

where we have used eqs. ( l24l) and ( |25l) . This can be recast as an upper bound on the neutrino 
Yukawa coupling y. Taking, e.g., the best fit values of the solar and atmospheric oscillation 
parameters given in Table d], we get: 

y < 0.035 (0.21) for NH with M x = 100 GeV (1000 GeV) and sin^ 3 = 0.1 , (30) 
y < 0.025 (0.15) for IH with M x = 100 GeV (1000 GeV) and sin 13 = 0.1 . (31) 

The constraints which follow from the current MEG upper bound on BR(/i — > e'y) will 
not be valid in the case of a cancellation between the different terms in one of the factors 
\Ue3 + i>ym2/m 3 Ue2\ 2 and \Ue2 + iUei\ 2 , I — e, [i, in the expressions (1271) and (1281) for BR(/i — > 
e'y). Employing the standard parametrisation of U, eqs. (12T1) - (T23]) . one can show that 
in the case of NH spectrum we can have |f/ e3 + i\/m 2 /m 3 U e 2\ = if [28] (see also |24j ) 
sin(5 + (a 2 i _ «3i)/2) = 1 and tan# 13 = (Ami/Am^ ) 1//4 sin6' 12 . Using the 3a allowed 
ranges of AttIq, Am\ and sin 2 # 12 given in Table (TJ we find that the second condition can 
be satisfied provided sin 2 #i 3 ^ 0.04, which lies outside the 3a range of allowed values of 
sin 2 #i 3 found in the Daya Bay experiment [11] (see eq. ([T|)). 

In the case of IH spectrum, the factor |C/ e2 + iU e \\ 2 can be rather small for sin(a 2 i/2) = 
— 1: \U e 2 + iU e i\ 2 = c 2 3 (l — sin26 l 12 ) = 0.0765, where we have used the best fit values of 
sin 2 9 12 = 0.306 and sin 2 # 13 = 0.0236. It is also possible to have a strong suppression of the 
factor \Ufj_2 + iU^ 2 EH]- Indeed, using the standard parametrisation of the matrix U it is 
not difficult to show that for fixed values of the angles #i 2 , 6* 23 and of the phases a 2 i and S, 
\U^2 + iU^ 2 has a minimum for 

c 23 cos 2#i 2 cos 5 sin ^f 1 — cos ^t 1 sin 5 , . 

s 23 1 + 2ci 2 si 2 sin ?f- 
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At the minimum we get: 



,, lON ( cos 5 cos + cos 29 12 sin 5 sin 221) 2 , 
mind^ + ^l 2 = 4 A 2 - 2J ~- (33) 

Notice that, from the equation above, the /1 — )■ e7 branching ratio is highly suppressed if 
the Dirac and Majorana phases take CP conserving values, mainly: 5 ~ and «2i — 7r. In 
this case, from eq. ( 132]) we get the lower bound sin(# 13 ) > 0.13, which is in agreement with 
the Daya Bay measurement reported in Tab. [TJ On the other hand, assuming CPV phases, 
we still may have min(|[/ M 2 + ^^i| 2 ) = 0, provided 6 , 12 and the Dirac and Majorana phases 
5 and «2i satisfy the following conditions: cos 5 cos(«2i/2) + cos 2#i 2 sin 5 sin(a 2 i/2) = 
and sgn(cos5cos ^f-) = — sgn(sin5sin ^-). Taking cos 5 > (cos 5 < 0) and using tan 5 = 
— tan(a 2 i/2)/ cos2# 12 in eq. (|32|) . we get the relation between si 3 , 5 and cos 29 12 , for which 
min(|f/ M 2 + ^i| 2 ) = 0: 



c 23 Vl + tan 2 <5 cos 26»i2 
sm6>i 3 = . (34) 

s 23 y 1 + cos 2 2# 12 tan 2 5 + 2c i2 s 12 sgn(cos 5) 

Using the 3a intervals of allowed values of sin 2 6*12 and sin 2 6*23 (found with the "new" reactor 
v e fluxes, see Tabled]) and allowing 5 to vary in the interval [0,27r], we find that the values 
of sin^ obtained using eq. ( 13~4"]) lie in the interval sin^ 13 ^ 0.11. As it follows from eq. ([[]), 
we have at 3<r: 0.10 ^ sin^ ^ 0.19. The values of 0.11 ^ sin^ ^ 0.19 correspond 
to < 6 ^ 0.7. These conclusions are illustrated in Fig. [1] For sin^i 3 and 5 lying in the 
indicated intervals we can have \U^2 + iU^i\ 2 = and thus a strong suppression of the fi — > 
decay rate. As we will see in subsections 2.2 and 2.3, in the model we are considering, the 
predicted fi — e conversion rate in a given nucleus and /1 — > 3e decay rate are also proportional 
to KRVy^RV) e i\ 2 , as like the fi — > e^y decay rate. This implies that in the case of the TeV 
scale type I see-saw mechanism and IH light neutrino mass spectrum, if, e.g., BR(/i — > cy) 
is strongly suppressed due to \U^2 + iU^\\ 2 = 0, the ji — e conversion and the fi — > 3e decay 
rates will also be strongly suppressed 0. The suppression under discussion cannot hold if, for 
instance, it is experimentally established that 5 is definitely bigger than 1.0. That would be 
the case if the existing indications [S] that cos 5 < receive unambiguous confirmation. 

The limits on the parameters K-RU)^! and |(i?V) e i|, implied by the electroweak precision 
data, eqs. ( JT8]) - ( 120]) . and the upper bound on BR(/i — > ej), eq. (|3]), are illustrated in Fig. 
[2] The results shown are obtained for the best fit values of sin #i 3 = 0.156 and of the other 
neutrino oscillation parameters given in Table [1] 



2.2 The /x — e Conversion in Nuclei 

We will discuss next the predictions of the TeV scale type I see-saw extension of the SM for 
the rate of the \i — e conversion in nuclei, as well as the experimental constraints that can 

5 Let us note that in the case of IH spectrum we are discussing actually one has |(i?F) M i| 2 oc 
\U^2 + «v /m i/ m 2^i| 2 ( see ec L- QUO; witn m 2 = \f\^- m W and m 1 = ^J\Am\\ - Am| = ^/|Amj[|(i - 
0.5Atoq/|Atoa|)- Therefore when |f/ M 2 + iU^i\ = we still have |£/ M 2 + i-\/ TO i/ m 2^i| 2 ^ 0. However, in 
this case |C/ m2 + i'\/ m i/ m 2^i| 2 — (Am|/(4|Am^|)) 2 |[/ A1 i| 2 1.7 x 1CT 5 , where we have used 6 = (which 
maximises |J7^i| ) and the best fit values of the other neutrino oscillation parameters. Thus, our conclusions 
about the suppression of BR(/i — > e^f), the /1 — e conversion and the (i — > 3e decay rates are still valid. 
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Figure 1: Values of siu $13, as a function of the phase 5, which yield a suppressed rate of the 
process \i — > ey. The values are obtained using eq. the 2a (3cr) intervals of allowed 

values of sin 2 $12 and sin 2 #23, yellow (green) points (found with the "new" reactor v e fluxes, 
see Tabled]) and allowing 6 to vary in the interval [0,27r]. The red and blue horizontal lines 
correspond to the 3a upper limit sin #13 = 0.191 and the best fit value sin 6*13 = 0.156. 



be imposed on this see-saw scenario by the current and prospective fi — e conversion data. 
In the type I see-saw scenario of interest, the fi — e conversion rate in a nucleus M is very 
well approximated by the expression [5"3] : 



CR(fiAf eAf) = 



"em e// 



capt 



2tt 4 



-m, 



r 



{Rvy^RvuY \c 



capt 



/A C I 



(35) 

In eq. fl35|) Z and are the proton number and the neutron number of the nucleus A/", 
respectively, F(— m 2 ) is the nuclear form factor at momentum transfer squared q 2 = — m 2 , 
being the muon mass, Z e ff is an effective atomic charge and r capt is the experimentally 
known total muon capture rate. The loop integral factor receives contributions from 
7— penguin, Z°— penguin and box type diagrams and reads [5 



m. 



fie — 



Z{AR {X) + H' (X))-(2Z + N) 



MX) 
sin 2 9 W 



(Z + 2N) 



Y (X) 
sin 2 9 



(36) 
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Figure 2: Correlation between |(_RV) e i| and |(i2V) M i| in the case of NH (upper panels) and 
IH (lower panels) light neutrino mass spectrum, for Mi = 100 (1000) GeV and, i) y — 0.0001 
(magenta points), ii) y = 0.001 (blue points), Hi) y = 0.01 (red points) and iv) y = 0.1 
(cyan points). The constraints from several LFV processes discussed in the text are shown. 



Here 6w is the weak mixing angle, sin 2 9w = 0.23, Z and iV are the proton and neutron 
numbers of the nucleus Af, X = Mf/M^, and 
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log y, 
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(37) 
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Figure 3: The fi — e conversion loop integration factor C^ e versus the see-saw mass scale M l5 
for three different nuclei: i) ^fTi (blue line), ii) ^Al, (green line), and Hi) 7g 7 Au (red line). 



In what follows we will present results for three nuclei which were used in the past, 
and are of interest for possible future \i — e conversion experiments: ^Ti, 13AI and 7g 7 Au. 
For these nuclei one has, respectively: i) Z e fj = 17.6; 11.62 ; 33.64, ii) F(q 2 = —mV) ~ 
0.54; 0.64; 0.20, and iii) T capt = 2.59; 0.69; 13.07 x lOSec" 1 [3]]. 

The dependence of the loop integration factor C Me on the see-saw mass scale M 1 for the 
three nuclei of interest is shown in Fig. The first feature to notice is that for &Ti, flAl 
and \f Au goes through zero at M 1 = 271; 285; 239 GeV, respectively. At M x = 271 GeV, 
\C„e\ for l 7 3 A\ and ^ 7 Au takes the values |C Me (Al)| ^ 8.6 and |C Me (Au)| ^ 158.0; at M x = 285 
GeV, we have |C Me (Ti)| = 16.2 and |C Me (Au)| ^ 234.5; and finally, at M x = 239 GeV, we 
find |C Me (Ti)| ^ 33.4 and |C Me (Al)] ^ 26.3. Thus, for the values of M x at which C Me (Ti) = 
or C M e(Al) = 0, |C Me (Au)| has a rather large value, while at the value at which C Me (Au) = 0, 
|C Me (Ti)| and |C Me (Al)| are not strongly suppressed. 

The interval where, e.g., |C Me | < 1 (10) in the case of fi — e conversion on Ti reads: 
270 (262) GeV £ |C Me | £ 272 (280) GeV. For a conversion on Al or Au these intervals are 
284 (269) GeV £ |C Me | £ 269 (301) GeV and 239.0 (237) GeV £ \C^\ £ 239.4 (241) GeV, 
respectively. 

For Mi lying outside these relatively narrow intervals centered on the values of M\ at 
which C^e = 0, \Cfj, e \ takes rather large values for each of the three nuclei. For M = 100 GeV, 
for instance, we find \C^ e (Ti)\ = 102, (^(Al)] = 61 and |C Me (Au)| ^ 367; for M = 1000 GeV 
one has: |C Me (Ti)| = 1936, |C^ e (Al)| = 1024 and \C^ e {Au)\ ^ 9110. 
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Figure 4: The ratio of the fi — e relative conversion rate and the branching ratio of the i) 
ix — > e~f decay (solid lines), ii) fx — > 3e decay (dashed lines), versus the type I see-saw mass 
scale Mi, for three different nuclei: ||Ti (blue lines), 13AI (green lines) and 7g 7 Au (red lines). 



It follows from these considerations that the \i — e conversion rate on, e.g., Ti, is strongly 
enhanced by the value of the loop integration factor \C^ e \ for M\ = (100 — 250) GeV and for 
Mi ^ 300 GeV. Similar results hold for the [i — e conversion rate on Al and Au nuclei. At 
Mi = 1000 GeV, for instance, the corresponding enhancement factors are ~ 10 6 for the Ti 
and Al nuclei [§, while for the Au nucleus it is ~ 10 7 . As a consequence, the rate of [i — e 
conversion has a remarkable sensitivity to the product |(i2V)* 1 (i?V) e i| of CC couplings of 
the heavy Majorana neutrinos to the electron and muon for most of the values of Mi from 
the interval of interest (100 — 1000) GeV, the only exception being the narrow intervals of 
values of Mi, quoted above, for which the conversion rate is suppressed. 

The best experimental upper bound on the conversion rate is [22]: CR(/iTi — > eTi) < 
4.3 x 10~ 12 . This bound implies a constraint on \(RV)* til (RV) e i\, which is shown in Fig. [2] 
for Mi = 100; 1000 GeV. It is quite remarkable that, as Fig. [2] shows, the constraint on 
the product of couplings \(RV)* fil (RV) e i\ implied by the best experimental upper limit on 
CR(/i Ti — y e Ti) is more stringent than the constraint following from the best experimental 
upper limit on BR(/i — > cy) although the two experimental upper limits are very similar 
quantitatively and the expression for CR(/iTi — > eTi) has an additional factor of a = 1/137 
with respect to the expression for BR(/i — > ej). 

6 Note that although IC^ej 2 increases monotonically with the increase of Mi ^ 300 GeV, the factor 
\(RV)* lll {RV) e i\ 2 in the expression for the conversion rate CR(fiAf —> ej\f) decreases as M^, so that 
CR(/iA/" — s- ej\f) has a finite limit when Mi tends to infinity. 
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Future experimental searches for fx — e conversion in ||Ti can reach the sensitivity of 
CR(/iTi — > eTi) ~ 10~ 18 [31]. Therefore, for values of Mi outside the narrow intervals 
quoted above for which the loop integration factor {C^l is strongly suppressed, an upper 
bound on the /i — e conversion ratio of O(10~ ls ) can be translated into the following stringent 
constraint on the heavy Majorana neutrino CC couplings to the muon and electron: 

\(RV)* fll (RV) el \ < 2.6 x 10" 8 (0.14 x 10" 8 ) for M 1 « 100 (1000) GeV . (41) 

As was noticed earlier, the two parameters of the type I see-saw model considered, the 
mass scale Mi and the Yukawa coupling y, can be determined, in principle, from data on 
BR(/i —> ej) (or BR(yU — > 3e)) and CR(/iTi — > eTi) if the two processes will be observed. 
Actually, the ratio of the rates of fi — e conversion in any given nucleus J\f, CR(yU M —> eAf), 
and of the [i — > e'j decay, depends only on the mass (scale) Mi and can be used, in principle, 
to determine the latter. In the case of fi — e conversion on titanium, for instance, we find: 

R ( -V—lS) = CR(^TweTi) _ 42 for M x « 100 (1000) GeV. (42) 

\fj, — > cy J BR(/i — > cy) 

The correlation between CR(fiAf —> eAf) and BR(/i — > ej) in the model considered 
is illustrated in Fig. HJ As Fig. |4] indicates, a measured value of R( ^ 5 would be 

compatible with two values of Mi. 

The type I see-saw mass scale Mi would be uniquely determined if \i — e conversion is 
observed in two different nuclei or if, e.g., the /i — > decay or fi — e conversion in a given 
nucleus is observed and it is experimentally established that R( ) ^ 10~ 3 . In the latter 
case Mi could be determined with a relatively high precision. Furthermore, if one restricts to 
RH neutrino masses in the range (50 — 1000) GeV, the mass scale Mi would also be uniquely 
determined if it is found that R( M ~ e ) > 20. 

We note also that the correlation between CR(/iA/" — > eAf) and BR(/i — > cy) in the 
type I see-saw model considered is qualitatively and quantitatively very different from the 
correlation in models where the fi — e conversion is dominated by the photon penguin dia- 
gram, e.g., the supersymmetric high-scale see-saw model which predicts approximately [56] 
CR(/iTi eTi) w 5 x 10~ 3 BR(/i ei). 



2.3 The v 3e decay 

The /i — > 3e decay branching ratio has been calculated in [51] in a supersymmetric scenario 
with lepton flavour violation. Adapting their results to the case we are interested in we get: 

4(T^ 9 9 

BR(/i^3e) = —f- ^RVy^RVU C, 3e (M 2 jM^) , (43) 

\CMv)\ 2 = 3Z\y) + 3Z(y)H(y) + H\y)(\og 1 ^-±f) 

\ m e 8 / 

+ ^?r Fe2 "Jr zye "Jr*' (44) 

A sin u\y sin t>w sin uw 

Here H{y) and Z{y) correspond to the dipole and charged radius type operators originating 
from the photon penguin diagram of the decay, while Y e (y) corresponds to the box diagram 
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Figure 5: The fi — > 3e decay rate factor | C M 3 e | 2 as a function of the see-saw mass scale Mi. 



contribution. We have: H(y) = (G(y) - G(0))/4, Y e (y) = Y (y) - \(RV) e i\ 2 S (y) « Y (y), 
where G(y) and Y (y) are given in eqs. ( )26l) and ( l40l) . respectively, and Z{y) = Co(y) + 
(D (y) - 2E {y)/3)/A with 



C (y) 
D (y) 



y-6 3y + 2 

y-i {y- 1) 2 

?/(18 + y - 7y 2 

E {y) = 



12(2/ - 1) £ 



y 2 (12-10y + y 2 
6(2/ -l) 4 



log 2/ • 



(45) 
(46) 



The dependence of the /i — > 3e decay rate factor (C^ael 2 on the type I see-saw mass 



scale Mi is shown in Fig. |5j At Mi 



fi3e | 

100 (1000) GeV we have: \C^ e \ 



1.1 (1639), i.e., 



|C M 3 e | increases by 3 orders of magnitude when Mi changes from 100 GeV to 1000 GeV. 
Using the quoted values of IC^I 2 we get the following constraint from the current limit on 
BR(/i 3e), eq. ®: 

\(RV)li(RV) e i\ < 2.1 x 10" 4 (5.3 x 10~ 6 ) for M x = 100 (1000) GeV. (47) 

Thus, for Mi = 100 GeV the constraint on \(RV)*^ ll (RV) e i\ obtained using the current 
experimental upper limit on BR(// — > 3e) is by a factor of 2.5 less stringent than that 
obtained from the current upper limit on BR(// — > ej) (see eq. ( 129]) ). while for Mi = 1000 
GeV it is by a factor of 6 more stringent. However, for the two values of Mi considered, the 
upper limit on \ (RV)*^i(RV) e i\ from the current experimental bound on the \i — e conversion 
rate, eq. §5§ is the most stringent. This is clearly seen in Fig. [2J It follows also from Fig. [2] 
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that an experiment sensitive to a /x — e conversion rate CR(/x Al — >■ e Al) « 10 16 , will probe 
smaller values of the product of couplings |(i?y)* 1 (i?l A ) e i| than an experiment sensitive to 
BR(/x -» 3e) = 10~ 15 . 

In Fig. H] we show the correlation between CR(/XjV — >■ c/V) and BR(/x — > 3e) in the 
TeV scale see-saw model considered. As it follows from Fig. [4] the observation of the 
/i — > 3e decay or of the \i — e conversion in a given nucleus, combined with data on the ratio 
CR(/xA/" —> eA/")/BR(/x — > 3e) would lead either to a unique determination of the type I 
see-saw scale Mi, or to two values, or else to a relatively narrow interval of values, of Mi 
compatible with the data. One can get the same type of information on the scale Mi from 
data on the ratio BR(/x — > 3e)/BR(/x — > e7), provided at least one of the two decays /x — > 
and n — > 3e is observed. 

It should be added finally that for Mi ^ 100 GeV we have: BR(/x — > 3e)/BR(/i — > 
ej) ^ 0.067. Thus, if it is experimentally established that BR(/x — > 3e)/BR(/x — > ej) is 
much smaller than the quoted lower bound, the model considered with Mi <^ 100 GeV will 
be ruled out. Such a result would be consistent also just with a see-saw scale Mi < 100 
GeV. 



3 TeV Scale Type II See-Saw Model 

We will consider in this section the type II see-saw [SZ] extension of the SM for the generation 
of the light neutrino masses. In its minimal formulation it includes one additional SU(2) l 
triplet Higgs field A, which has weak hypercharge Yy/ = 2: 

A ~ { A° -A+/V2 ) ■ (48) 

The Lagrangian of the type II see-saw scenario, which is sometimes called also the "Higgs 
Triplet Model" (HTM), reads 

Cesaw = -M 2 A Tr (At A) - (h u > WiL *T2 A # £ + M H T tr 2 A^H + h.c.) , (49) 

where (^l) t = iyjt ^l), ip c \l = ( — v 1lP~ X — ^iP _1 ), and H are, respectively, the 
SM lepton and Higgs doublets, C being the charge conjugation matrix, and /ia is a real 
parameter characterising the soft explicit breaking of the total lepton charge conservation. 
We are interested in the low energy see-saw scenario, where the new physics scale Ma 
associated with the mass of A takes values 100 GeV < Ma < 1 TeV, which, in principle, 
can be probed by LHC [59] . 

The flavour structure of the Yukawa coupling matrix h and the size of the lepton charge 
soft breaking parameter /xa are related to the light neutrino mass matrix m„, which is 
generated when the neutral component of A develops a "small" vev v& oc /xa • Indeed, 
setting A = va and H T = (0 v) T with v ~ 174 GeV, from Lagrangian ( 14*91) one obtains: 

(m v ) u , = mw - 2 hw va ■ (50) 



7 We do not give here, for simplicity, all the quadratic and quartic terms present in the scalar potential 
(see, e.g., [55]). 
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The matrix of Yukawa couplings hw is directly related to the PMNS neutrino mixing matrix 
[TpMNS = U, which is unitary in this case: 



hw = 2^ ( U * diag(mi, m 2 ,m 3 ) U ] ) u , . (51) 

An upper limit on v A can be obtained from considering its effect on the parameter p = 
M^/Ml cos 2 6 W . In the SM, p = 1 at tree-level, while in the HTM one has 

1 + 2x 2 

p=l + dp= 1 + 4x2 , x = v A /v. (52) 



The measurement p ~ 1 leads to the bound v A /v < 0.03, or v A < 5 GeV (see, e.g., [60]). 

As we will see, the amplitudes of the LFV processes p, — > ej, p — > 3e and p + Af —> e + Af 
in the model under discussion are proportional, to leading order, to a product of 2 elements of 
the Yukawa coupling matrix h. This implies that in order for the rates of the indicated LFV 
processes to be close to the existing upper limits and within the sensitivity of the ongoing 
MEG and the planned future experiments for M A ~ (100 — 1000) GeV, the Higgs triplet 
vacuum expectation value v A must be relatively small, roughly v A ~ (1 — 100) eV. In the 
case of M A ~ v = 174 GeV we have v A = p, A , while if M\ » v 2 , then v A = p, A v 2 /(2M A ) 
(see, e.g., [SSI ED]) with v 2 j(2M\) ^ 0.015 for M A = 1000 GeV. Thus, in both cases a 
relatively small value of v A implies that p, A has also to be small. A nonzero but relatively 
small value of p, A can be generated, e.g., at higher orders in perturbation theory [61] or in 
the context of theories with extra dimensions (see, e.g., [62]). 

The physical singly-charged Higgs scalar field practically coincides with the triplet scalar 
field A + , the admixture of the doublet charged scalar field being suppressed by the factor 
v A /v. The singly- and doubly- charged Higgs scalars A + and A ++ have, in general, different 
masses [SU [63]: m A + ^ m A ++. Both situations m A + > m A ++ and m A + < m A ++ are possible. 
In some cases, for simplicity, we will present numerical results for m A + = m A ++ = M A , but 
one must keep in mind that m A + and m A ++ can have different values. 

In the mass eigenstate basis, the effective charged lepton flavour changing operators arise 
at one-loop order from the exchange of the singly- and doubly-charged physical Higgs scalar 
fields. The corresponding effective low energy LFV Lagrangian, which contributes to the 
pi — e transition processes we are interested in, can be written in the form: 

C efS = ~4^f (m ll A R ea^P R p J F Pa + h.c.) 



V2 
V2 



A L {-ml)erPLV QQQlaQ + h.c.) , (53) 

Q=u,d 



where e is the proton charge, and q u = 2/3 and qd = —1/3 are the electric charges of the up 
and down quarks (in units of the proton charge). We obtain for the form factors Ar^: 

A R = - L { 1olr (At- + > (54) 



y/2G F 48vr 2 \8m 2 A+ m A+H 



Mq2) ~ V~2G F k t 2 (l2< + + m A++ ; (m^' ' ' ' (55) 
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mi being the mass of the charged lepton I, I = e,/j,,T. In the limit where the transition 
is dominated by the exchange of a virtual doubly charged scalar A ++ , these expressions 
reduce to those obtained in [SU [65]; to the best of our knowledge the expression of A^q 2 ) 
for the general case is a new result. The term with the form factor Ar in eq. (1531) generates 
the /jL — > e7 decay amplitude. It corresponds to the contribution of the one loop diagrams 
with virtual neutrino and A + [66J and with virtual charged lepton and A ++ [M] (see also 
[68]). The second term involving the form factor A^, together with A^, generates the fi — e 
conversion amplitude. The loop function f(r,si) is well known [65] : 



f(r,si) = — + log(sj) + 1 M/1 + — log^ -, ; . ■ (56) 



r \ r J V r y/r — \fr + 4s/ 

Notice that in the limit in which the charged lepton masses m/ are much smaller than 
the doubly-charged scalar mass m&++, one has f(r,si) ~ log(r) = \og(m 2 l /m 2 ± ++ ). For 
m&++ = (100 — 1000) GeV, this is an excellent approximation for f(r,s e ), but cannot 
be used for f(r, s^) and f(r,s T ). The ratios fir^e)/}^^^) and f(r,s e )/f(r,s T ) change 
relatively little when m^++ increases from 100 GeV to 1000 GeV, and at m&++ = 100 (1000) 
GeV take the values: /(r, s e )/f(r, s^) = 1.2 (1.1) and f(r,s e )/f(r,s T ) ^ 2.1 (1.7). More 
generally, f(r, si), I = e, /i, r, are monotonically (slowly) decreasing functions of m&++ [^: for 
m A++ = 100 (1000) GeV we have, e.g., /(r,s e ) = -13.7 (-18.3). 



3.1 The fi — > Decay 



Using eqs. ( 153"]) and ( 15"4"|) we get for the \i — > cy decay branching ratio in the case under 
discussion 1661 [581: 



BR(/i erf) = 384 vr 2 (4tt a cm ) \A R \ 



e/i 



192 7T 



m, 



A+- 



(57) 



For m&+ — m&++ = Ma, the upper limit on BR(/i — > ej) reported by the MEG experiment, 
eq. ([3]), implies the following upper bound on \(h^h) efJi \: 



(rfh) 



e/i 



< 5.8 x 10" 



Ma 



100 GeV J 



V. 



(55 



One can use this upper bound, in particular, to obtain a lower bound on the vacuum expec- 
tation value of A , v& 0. Indeed, from eq. (15T1) it is not difficult to get: 



(h)h) 



eft, 



U e2 UlAm 2 21 



U e3 Ul Am 2 31 



(59) 

where we have used the unitarity of U. The above expression for \(h^h) e ^\ is exact. It 
follows from eq. (1591) that the prediction for \(h^h) efl \, and thus for BR(/i — > ej), depends, 
in general, on the Dirac CPV phase 5 of the standard parametrisation of the PMNS matrix 
U (see eq. ( I2~il ). For the best fit values of sin 2 #13 = 0.0236 and of the other neutrino 



8 Note that we have f(r, si) < 0, I = e, fi, r. 
9 This was noticed also in 1571. 
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Figure 6: The dependence of \{h)h) ei± \ on sin 6^3 for v A = 9.5 eV and 5 = (solid lines) 
and 5 = ir/2 (dashed line). The other neutrino oscillation parameters are set to their best 
fit values given in Table [TJ The vertical line corresponds to the current 3a allowed minimal 
value of sin 9 13 (see eq. §B)). 



oscillation parameters listed in Table [H the term oc Am^ in eq. (159]) is approximately a 
factor of 10 smaller than the term oc Am^. In this case BR(/i — > cy) exhibits a relatively 
weak dependence on the type of the neutrino mass spectrum and on the Dirac phase S. 
Neglecting the term oc Am^, we obtain from (J5SJ) and (15 90 : 



v A >2.1x 10 2 |s 13 s 2 3 Am^ 



100GeV\ 
M A ) 



3.0 eV 



100GeV\ 



Ma 



(60) 



For the 3a allowed ranges of values of sin 2 2# 13 given in eq. ([T]) and of the other neutrino 
oscillation parameters quoted in Table [U the absolute lower bound on v A corresponds approx- 
imately to v A > 1.5 eV (100 GeV)/MA and is reached in the case of Am^ > (Am^ < 0) 
for 5 = vr (0). 

We note further that if 5 = n/2 (37r/2), the term oc Am^ in the expression for \{h)h) etl \ 
(and thus for BR(/i — > ey)) always plays a subdominant role as long as the other neutrino 
oscillation parameters lie in their currently allowed 3a ranges. Therefore in this case the 
dependence of BR(/i — > ej) on the type of neutrino mass spectrum is negligible. The specific 
features of the predictions for \{h)h) eiJ \ discussed above are illustrated in Fig. [6j 

Exploiting the fact that vW{h)h) e)J \ is known with a rather good precision, we can write: 



BR(/i -»■ 67) 2.7 x 10 



-10 



leV 



100 GeVV 
M A ) 



(61) 
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where we have used eq. (|57|) and the best fit values of the neutrino oscillation parameters. 
It follows from eq. (I61~j) that for the values of va and Ma (or f«A+ and/or ?tia++) of in- 
terest, BR(/i — > e7) can have a value within the projected sensitivity of the ongoing MEG 
experiment. 



3.2 The 11 ->> 3e decay 

In the TeV scale type II see-saw scenario, the fi — > 3e decay amplitude is generated at the 
tree level by the diagram with exchange of a virtual doubly-charged Higgs boson A ++ . The 
branching ratio can be easily calculated (see, e.g. 



BR /i -> 3e = —l^J^J^L = ' f| f , 62 

G| m^ ++ G 2 F m A++ 16 v% 

where we have used eq. (1501 . From the present limit BR(/i — >■ 3e) < 10~ 12 , one can obtain 
the following constraint on \{h + ) ee {h) ^ e \: 

2 



\(h%(hU < 1.2 xlO- 7 (y^^) • (63) 



In the model under discussion, BR(/i — > 3e) depends on the factor |m* e ra Me |, which 
involves the product of two elements of the neutrino Majorana mass matrix, on the neutrino 
mass spectrum and on the Majorana and Dirac CPV phases in the PMNS matrix U. For 
the values of iba+ and mA++ in the range of ~ (100 — 1000) GeV and of i>a <§C 1 MeV of 
interest, m ee practically coincides with the effective Majorana mass in neutrinoless double 
beta ((PP)ov-) decay (see, e.g., pi ESI ED]), (m): 



3=1 



(m) | . (64) 



Depending on the type of neutrino mass spectrum, the value of the lightest neutrino mass and 
on the values of the CPV Majorana and Dirac phases in the PMNS matrix, \m ee \ can take 
any value between and tuq, where m = mj = m 2 = m 3 is the value of the neutrino masses 
in the case of quasi-degenerate (QD) spectrum, m ^ 0.1 eV (see, e.g., [OH])- It follows from 
the searches for the (/3/3)o I /-decay that \m ee \ ^ mo ^ 1 eV, while the cosmological constraints 
on the sum of the neutrino masses imply mo ^ 0.3 eV (see, e.g., [1]). As is well known, 
the (/?/?) ovdecay is claimed to have been observed in [7T1 E2], with the reported half-life 
corresponding to [72] \m ee \ = 0.32 ± 0.03 eV. This claim will be tested in a new generation 
of (/3/3)o^-decay experiments which either are already taking data or are in preparation at 
present (see, e.g., [T| 173]). 

In the case of NH light neutrino mass spectrum with ni\ <^ 10~ 4 eV, \m ee \ lies in the 
interval 3.6 x 10" 4 eV £ \m ee \ ^ 5.2 x 10~ 3 eV. This interval was obtained by taking into 
account the 3a allowed ranges of values of sin 2 # 13 (eq. Q), sin 2 # 12 , sin 2 623 and Anig and 
Am 2 A . For the best fit values (b.f.v.) of the latter we geto 1.45 x 10~ 3 eV ^ \m ee \ ^ 3.75 x 



10 The numerical values quoted further in this subsection are obtained for the indicated best fit values of 
the neutrino oscillation parameters, unless otherwise stated. 
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Figure 7: The dependence of \m* ee m^ e \ on the lightest neutrino mass mi in the case of 
neutrino mass spectrum with normal ordering (Am\ > 0), for four sets of values of the 
Dirac and the two Majorana CPV phases, [5,021,0:31]. The depicted curves correspond to 
the best fit values of sin #13 (eq. ([I])) and of the other neutrino oscillation parameters given 
in Table HJ The scattered points are obtained by varying the neutrino oscillation parameters 
within their corresponding 3cr intervals and giving random values to the CPV Dirac and 
Majorana phases. 



10 -3 eV. The minimal and the maximal values correspond to the combination of the CPV 
phases (021 — 031 + 25) = it and 0, respectively. However, for vii\ ^ 10~ 4 eV, one can have 
\m ee \ = for specific values of mi if the CPV phases 021 and 031 — 25 possess the CP 
conserving values 021 = vr and (031 — 25) = 0,tt (see, e.g., [H]): for the [ir,0] combination 
this occurs at m\ = 2.3 x 10~ 3 eV, while in the case of the [vr, 7r] one we have \m ee \ = 
at mi = 6.5 x 10~ 3 eV. If the light neutrino mass spectrum is with inverted ordering 
(Am\ = Amg 2 < 0, 7713 < mi < m2) or of inverted hierarchical (IH) type (m3 <C mi < m2), 
we have [71] |m ee | ^ \J \ Am\ | + m| cos 2Q\ 2 ^ 1-27 x 10~ 2 eV, while in the case of QD 
spectrum, \m ee \ ^ mo cos 2#i2 ^ 2.8 x 10~ 2 eV, where we used the 3a minimal allowed values 
of |Am^J and cos2#i2. 

We consider next briefly the dependence of the neutrino mass matrix element Im^l on 
the type of the neutrino mass spectrum and on the CPV Majorana and Dirac phases. In the 
case of NH spectrum with mi = 0, the maximal value of \m^ e \ is obtained for 031 — «2i = 5, 
5 = 7r, and reads: maxdm^l) = 8.1 x 10~ 3 eV. We get (m^l = for o 2 i = tt, 5 = (n) and 
«3i = (tt). As can be shown, for each of these two sets of values of the CPV phases, the zero 
takes place at essentially the same value of mi = 4.3 x 10 -3 eV (Fig. [7]). If the neutrino mass 
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Figure 8: The same as in Fig. [7] in the case of a light neutrino mass spectrum with inverted 
ordering (Am\ < 0) (see text for details). 



spectrum is of the IH type with negligible m 3 = 0, the maximal value of {m^l corresponds 
to 5 = and a 2 i = ^ and is given by maxdm^l) = \J \ Am\ | ci 3 (c 23 sin 26 12 + s 23 Si 3 cos 26*12) • 
The element \m^ e \ is strongly suppressed, i.e., we have |m Me | <C maxdm^l), for 5 = ir/2 and 
a value of the Majorana phase a 2 i which is determined by the equation: 

C23 C12 S12 sin a 21 = (c\ 2 + s 2 2 cos a 2 i) s 23 S13 • (65) 

For the b.f.v. of the neutrino mixing angles this equation is satisfied for a 2 \ — 0.283. 

The properties of \m ee \ and \m^ e \ described above allow us to understand most of the 
specific features of the dependence of the quantity \m* ee m^ e \ of interest on the the neutrino 
mass spectrum and the leptonic CPV phases. For NH spectrum and negligible mi = 0, the 
maximum of the latter is obtained for a 3 i — a 2 \ = 5 = and is given by: 

max( I m* e m Me |) = | (m 2 s\ 2 c? 3 + m 3 s 2 13 ) c 13 (m 2 s 12 (c 12 c 23 - s 12 s 23 s 13 ) + m 3 s 23 s 13 )\ , 

(66) 

with m 2 = ^ArriQ and m 3 = */ Am\. Using the b.f.v. of the neutrino oscillation parameters 
we get max(|m* e m Me |) = 2.9 x 1(T 5 eV 2 (see Fig. EJ). This implies BR(/i -> 3e) ^ 6 x 
10~ 9 (1 eV/vA) 4 (100 GeV/mA++) 4 . In the case of NH spectrum and non-negligible m x we 
have \m* ee m^ e \ = for the values of the CPV phases and m\ discussed above, for which either 
\m ee \ = or jm^el = 0. The scattered points in Fig. [7] correspond to the possible values 
the quantity |m* e m Me | can assume when varying the neutrino oscillation parameters within 
their corresponding 3a intervals and giving random values to the CPV Dirac and Majorana 
phases from the interval [0,27r]. 
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The maximum of \m ee m^ e \ for the IH spectrum with a negligible m 3 is reached for 5 = 
and CC21 = 7T, and reads: 

max( I m* ee m^ e |) ^ |Am|| c 3 3 Q c 23 sin46»i 2 + s 23 S13 cos 2 26»i 2 ^ . (67) 

Numerically this gives max(\m ee m^ e \) = 6.1 x 10~ 4 eV 2 (Fig. [8]). For BR(/i — > 3e) we thus 
obtain: BR(/i 3e) £ 2.4 x 10~ 6 (1 eV/w A ) 4 (100 GeV/m A ++) 4 . One can have \m ee m^ e \ <C 
max(|m ee m Me |) in the case of IH spectrum with m 3 = for, e.g., 5 = it/2 and a 2 i = 0.283, 
for which {m^l has a minimum. For the indicated values of the phases we find: |m ee m Me | = 
1.2 x 10~ 6 eV 2 (see Fig. |H]). Similarly to the case of a neutrino mass spectrum with normal 
ordering discussed above, we show in Fig. [8] the range of values the LFV term Im^m^e! can 
assume (scattered points). 

Finally, in the case of QD spectrum, |m ee m Me | will be relatively strongly suppressed 
with respect to its possible maximal value for this spectrum (i.e., we will have \m ee m iie \ <C 
max^m^m^!)) if, e.g., the Majorana and Dirac phases are zero, thus conserving the CP 
symmetry: a 21 = a 31 = 5 = 0. Then one has: \Tn ee m^ e \ = | Am^|si 3 s 23 Ci 3 /2 = 1.2 x 
10 -4 eV 2 . Note that this value is still larger than the maximal value of \m ee m^ e \ for the 
NH neutrino mass spectrum with a negligible mi (see Fig. [7]). The maximum of jr/wm^l 
takes place for another set of CP conserving values of the Majorana and Dirac phases: 
«2i = « 3 i = 71 an d 5 = 0. At the maximum we have: 

max(|m ee r?v|) ^ m 2 (c 3 3 cos26>i 2 - s 2 3 ) c 13 (c 23 sin26> 12 + 2c 2 2 s 23 s 13 ) , m ^ 0.1 eV. 

(68) 

For the b.f.v. of the neutrino mixing angles we get m^{\m* ee m^ e \) = 0.3m 2 ,. For mo ^ 0.3 
eV this implies max(|m* e m Me |) ^ 2.7 x 10~ 2 eV 2 , leading to an upper bound on BR(/i — > 3e), 
which is by a factor approximately of 4.1 x 10 3 larger than in the case of IH spectrum. 
The features of \m ee m^ e \ discussed above are illustrated in Figs. [7] and [81 

It should be clear from the preceding discussion that in the case of the type II see- 
saw model considered, the value of the quantity |(^) ee (^)^e| 2 oc |m* e m Me | 2 , and thus the 
prediction for BR(/i — > 3e), depends very strongly on the type of neutrino mass spectrum. 
For a given spectrum, it exhibits also a very strong dependence on the values of the Majorana 
and Dirac CPV phases a 2 i, a 3 i and 5, as well as on the value of the lightest neutrino mass, 
min(m,'). As a consequence, the prediction for BR(yU — > 3e) for given v& and m A ++ can 
vary by a few to several orders of magnitude when one varies the values of min(mj) and 
of the CPV phases. Nevertheless, for all possible types of neutrino mass spectrum - NH, 
IH, QD, etc., there are relatively large regions of the parameter space of the model where 
BR(/i — y 3e) has a value within the sensitivity of the planned experimental searches for the 
fi — > 3e decay [33]. The region of interest for the NH spectrum is considerably smaller than 
those for the IH and QD spectra. In the case NO spectrum (Am\ > 0), BR(/i — > 3e) can 
be strongly suppressed for certain values of the lightest neutrino mass mi from the interval 
~ (2 x 10~ 3 — 10~ 2 ) eV (Fig. [7J). For the 10 spectrum (Am\ < 0), a similar suppression 
can take place for m 3 10~ 2 eV (Fig. [8]). In the cases when |m* e m Me | 2 is very strongly 
suppressed, the one-loop corrections to the fi — > 3e decay amplitude should be taken into 
account since they might give a larger contribution than that of the tree level diagram we 
are considering. The analysis of this case, however, is beyond the scope of the present 
investigation. 
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48 Ti 

22 11 


0.0864 


0.0396 


0.0468 


2.590 




0.0362 


0.0161 


0.0173 


0.7054 


7 fAu 


0.189 


0.0974 


0.146 


13.07 



Table 2: Nuclear parameters related to \i — e conversion in ffjTi, ^Al and 7 g 7 Au. The 
numerical values of the overlap integrals D, and are taken from [75] . 



3.3 The fi — e Conversion in Nuclei 



Consider next the fi—e conversion in a generic nucleus A/". We parametrise the corresponding 
conversion rate following the effective field theory approach developed in [75]. Taking into 
account the interaction Lagrangian ( 153]) . we get in the type II see-saw scenario 



CR(/iA/" -)• eJV) 
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(69) 



The parameters D and represent overlap integrals of the muon and electron wave func- 
tions and are related to the effective dipole and vector type operators in the interaction 
Lagrangian, respectively (see, e.g. [75]). 

In the case of a light nucleus, i.e. for Z < 30, one has with a good approximation 
D ~ 8 \fATva cra V^ p \ with the vector type overlap integral of the proton, V^ p \ given by 

1 



m 



5/2 3/2 y2 7 l/2 p/ 



-m 



(70) 



where F(q 2 ) is the form factor of the nucleus. The parameters D m^ 5 ^ 2 , m At 5//2 and r capt 
for 22 Ti, 13AI and 7 9 7 Au are given in Table [2j Using the result for D quoted above, eqs. (154]) . 
(155]) and ( 170]) . the conversion rate (169]) can be written as 
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(71) 



In contrast to previous studies, which assume that the /1 — e conversion is dominated by the 
A ++ exchange (65] , we will consider in this work a scenario where both scalars contribute to 
the transition amplitude. Thus, assuming m A + = m A ++ = Ma, we have CR(/iA/" —> eAf) oc 
|C^ 7) | 2 , where 



C (ii) 

w fie 



1 

4^1 



29 , t v 
— (m 1 m 



mjj(r,s0m^ , (72) 

l=e,/i,T 

and we have used eq. (1501 . The upper limit on the fi — e conversion rate in Ti, eq. ([5]), leads 
to the following upper limit on \C\!P\: 

X. 
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(73) 
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Figure 9: The dependence of 4v^|C^ e | (given in eV 2 ) on the lightest neutrino mass m\ in 
the case of neutrino mass spectrum with normal ordering (Am| > 0), for two sets of values 
of the Dirac and the two Majorana CPV phases, [5, an, 0:31] and Ma = 200 (1000) GeV, 
plain (dashed) curves. The figure is obtained for the best fit values of sin #13 (eq. (JTJ) and 
of the other neutrino oscillation parameters given in Tabled] (see text for details). 



In obtaining it we have used the values of r capt , Z e ff, Z and F{— m 2 ) for Ti given in 
subsection 2.2. An experiment sensitive to CR(/iTi — > eTi) « 10~ 18 [31] will be able to 
probe values of \C$P\ £ 5.8 x 10~ 8 (M A /(100 GeV)) 2 . 

The /i — e conversion rate in a given nucleus depends through the quantity C^J\ on the 
type of neutrino mass spectrum and the Majorana and Dirac CPV phases in the PMNS 
matrix. Using the b.f.v. of the the neutrino oscillation parameters and performing a scan 
over the values of the CPV phases and the lightest neutrino mass, which in the cases of 
NO (Am} > 0) and IO (Am} < 0) spectra was varied in the intervals (10~ 4 — 1) eV and 
(10~ 5 — 1) eV, respectively, we have identified the possible ranges of values of Av\\C\!P\. 
The latter are shown in Figs. (Q and f fTOj) . 

For M A = 200 (1000) GeV and NH spectrum with negligible m 1 (m 1 < 10" 3 eV), the 
maximal value of Av\\C\fP\ occurs for [5, «2i,«3i] = [0,0,0] and at the maximum we have 
Av 2 A \C ( J e I] \ = 2.9 (3.8) x 10~ 3 eV 2 . For values of the CPV phases [5,a 21 ,a 3 i] = [0, tt, 0] and 
M A = 200 GeV, 4t> A |c£ 7) | goes through zero at m 1 = 2 x 10~ 2 eV (Fig. EJ). In the case of 
a larger charged scalar mass, i.e. Ma = 1000 GeV, such cancellation occurs at a different 
value of the lightest neutrino mass, mainly m\ = 0.025 eV. 

The maximum of Av\\C^!P\ in the case of IH spectrum with negligible m 3 , occurs for 
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Figure 10: The same as in Fig. [9]in the case of a light neutrino mass spectrum with inverted 
ordering. 



maximal CPV phases: [8, 021 , a^x] — [vr/2, 37r/2, 0]. At the maximum in this case one has 
4w A |c£f ) | ^ 6 (7) x 10" 3 eV 2 for M A = 200 (1000) GeV. As Fig. M shows, for other sets 
of values of the CPV phases, Av\\C^P\ can be much smaller. Taking again CP conserving 
phases, e.g. [iv,iv, 0], one can get a strong suppression of the branching ratio for m 3 = 
7.2 (15) x 10" 3 eV and M A = 200 (1000) GeV. Allowing sin6>i 3 to take values other than the 
best fit one, we find that Av\\C\!P\ can even go through zero at, e.g., [5, «2i> 031] = [tt, 7r, tt/2] 
for sin 6*i3 = 0.137, which lies within the 2a allowed region. In Fig. [10] we report other 
examples in which the CPV phases in the PMNS matrix take different sets of CP violating 
values and the quantity Av\\C^P\ (and the conversion rate) can vary by several orders of 
magnitude for specific values of the lightest neutrino mass 7713 and the see-saw mass scale 
M A . 

If the neutrino mass spectrum is quasi-degenerate, m 12 .3 = m ^0.1 eV, we have for 
m ^ 0.3 eV: 2.8 x 10" 3 eV 2 ^ Av A \C^\ ^ 0.4 eV 2 . The minimal value corresponds to 
Am^ > (NO spectrum) and [5, 021,031] = [tt, 0, 0]; for e.g. [5,021,031] = [0,0,0] and 
M A = 200 GeV we get in the QD region 4u!|Cf e 7) | 3.3 x 10~ 3 eV 2 (Fig. EJ. 

Finally, the scattered points in Figs. M and [10] are obtained by varying all the neutrino 
oscillation parameters within the corresponding 3a intervals and allowing for arbitrary values 
of the Dirac and Majorana phases in the interval [0,27r]. 

We remark that the previous estimates, as well as Figs M and [TU] were realized under 
the assumption that the singly- and doubly-charged scalars have masses of the same order, 
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i.e. — wa++ = Ma- The case in which the dominant contribution to the conversion 
amplitude is provided by the exchange of A ++ , i.e. for 3> m^++ > 100 GeV, shows 
similar features: the upper limits of the conversion ratio in the cases of NO and 10 spectra 
are unchanged and a strong suppression can occur for specific values of the CPV phases 
and min(mj). Taking, instead, the opposite limit m^++ ^> m^+, with = (100 — 1000) 
GeV, the dominant contribution to the fi — e conversion amplitude is given by the exchange 
of the singly-charged scalar, therefore we have: \C\iP\ oc \{h)h) eiJ \. As it was pointed out 
in subsection 13. 1[ \(h^h) efJ ,\ shows a relative weak dependance on the type of neutrino mass 
spectrum and on the CPV phases in the PMNS matrix. Moreover, no suppression of the 
conversion amplitude occurs if sin(#i3) is taken within the current 3cr experimental bound 
(see Fig. [6]). In this case, from the best experimental upper bound on the conversion rate in 
Ti, CR(/xTi -> eTi) < 4.3 x 10" 12 , we get the constraint: 

\{h)h) m \ < 6 x 1(T 4 ( m * + V , (74) 
lv ; Ml V 100 GeV/ ' V ; 

which provides a weaker bound with respect to that obtained from the /x — > e'j decay (see 

eq. (1581) ). A \x — e conversion experiment sensitive to i.e. CR(/iTi — > eTi) ps 10~ 18 , can 

probe values of \{h)h) e ^\ which are by a factor 2 x 10 3 smaller and could set the limit: 

'< fttft W < 3 * 10 " (l^v) 2 ' < 75 > 
4 TeV Scale Type III See-Saw Model 

We turn in this section to the study of lepton flavour violating processes in type III see- 
saw [76J extensions of the SM. In the scenarios under discussion, the SM particle content is 
enlarged by adding S77(2)£-triplets of fermions, FjR = {Fj R , Fj R , Fj R ) , j > 2, possessing 
zero weak hypercharge and a mass at the electroweak scale: ps (100 — 1000) GeV. 
The corresponding interaction and mass terms in the see-saw Lagrangian read: 

Cesaw = "A« tflL-rH- FjR - \ {Mr\. Jg ■ FjR + h.c. , (76) 

where r = (-T 1 , r 2 , r 3 ), r a being the usual SU{2)l generators in the fundamental represen- 
tation. 

It is convenient in the following discussion to work with the charge eigenstates Fp R = 
(Fj R =F iFjR) an d FjR = FjR- Then, the physical states in the above Lagrangian correspond 
to electrically charged Dirac and neutral Majorana fermions, which are denoted as Ej and 
Nj, respectively: o 

Ej ee F7 R + Fff Nj ee F% + Ff R . (77) 

In the basis in which the charged lepton mass matrix is diagonal, the CC and NC weak 
interaction Lagrangian of the light Majorana neutrino mass eigenstates Xj read: 

C v c c = - -£= l la ((1 - r 1 )U) H XiL W a + h.c. , (78) 
Cnc = -^-Wla(uHl + 2 V )U).. Xj LZ a . (79) 

In the following we will denote as Ej and Nj the mass eigenstates obtained from the diagonalization of 
the full charged and neutral lepton mass matrices. 
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Similarly to the type I see-saw scenario discussed earlier, the heavy Majorana mass eigen- 
states N k might acquire a sizable coupling to the weak gauge bosons through the mixing 
with the light Majorana neutrinos: 



C-cc = ^y=£ja(RV) ik (l-j 5 )N k W a + h.c, (80) 
C N NC = -^u ela (RV) e ,(l- l5 )N,Z a + h.c. (81) 

In the expressions given above, the non-unitary part of the neutrino mixing matrix, i.e. 
the matrix rj, and the matrix R are defined as in the type I see-saw scenario discussed in 
Section [2] (see eq. f)13p ). while V in this case diagonalizes the symmetric mass matrix M# in 
eq. CED: M R V*diag(Mi, M 2 , . . .)Vt 

The neutrino Yukawa couplings Xej can be partially constrained by low-energy neutrino 
oscillation data and electroweak precision observable (see, e.g. [771 [78]). Notice that, unlike 
the type I see-saw extension of the Standard Model, now we have flavour changing neutral 
currents (FCNCs) in the charged lepton sector. The latter are described by the interaction 
Lagrangian: 

4c = tt- (hlail-^u,^ -2sll la i) Z a . (82) 

Finally, the interactions of the new heavy charged leptons, Ej, with the weak gauge bosons 
at leading order in the mixing angle between the heavy and the light mass eigenstates read: 

L E CC = -gE jlot N j W a + gE jla (RV) ij ugW a + h.c, (83) 
C E NC = gc^Ej^EjZ" - (Z 7q (W)^-(1-7 5 )£,-Z q + h.c.) . (84) 

4.1 The fi — > e7 Decay 

Charged lepton radiative decays receive additional contributions with respect to the scenario 
with singlet RH neutrinos, due to the presence of new lepton flavour violating interactions 
in the low energy effective Lagrangian (see eqs. (1821) and (1841) ). Following the computation 
reported in [78], we have for the \i — > e^ decay branching ratio in the present scenario: 

BR(Me7) = ^\T\\ (85) 

OZ7T 

where the amplitude T is given by 

T = -2 + cj r)» e + (RV)ek(RV); k [A(x k ) + B(y k ) + C(z k ) ] , (86) 



12 Flavour changing couplings between the charged leptons and the SM Higgs boson H arise as well in the 
TeV-scale type III see-saw scenarios [75] which enter at one-loop in the lepton flavour violating processes 
(see next subsection). 
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with x k = (M k /M w ) 2 , y k = (M k /M z ) 2 , z k = (M k /M H ) 2 and C ~ -6.56. The loop functions 
A(x k ), B(y k ) and C(z k ) read [78]: 

-30 + 153x - 198a; 2 + 75x 3 + 18(4-3x)x 2 logx , , 

^ = 3(^1? " " ' (8?) 

33 - 18 y - 45y 2 + 30y 3 + 18 (4 - 3y) y logy 

Biv) = • ^ 

-7 + 12 z + 3z 2 - 8z 3 + 6(3,z-2)zlogz 
[Z) ~ 3(^-l) 4 • 1 j 

In the simple scenario of degenerate fermion triplets with an overall mass scale M we obtain 
taking M H = 125 GeV: 

T/^ e = 11.6(5.2), for M = 100 (1000) GeV . (90) 

For M = 100 (1000) GeV, the current best upper limit on the \i — > cy decay branching ratio 
obtained in the MEG experiment, eq. (J3]), implies the bound: 

l^el < 9 (20) x 10~ 6 , for M = 100 (1000) GeV. (91) 

If no positive signal will be observed by the MEG experiment, that is if it results that 
BR(/i — > &y) < 10~ 13 , the following upper limit on the non-unitarity lepton flavour violating 
coupling | rj^e | can be set: 

|^ e | < 2 (4) x 10" 6 , for M = 100 (1000) GeV. (92) 
4.2 The ji — > 3 e and /i — e Conversion in Nuclei 

The effective \x — e — Z effective coupling in the Lagrangian (182]) provides the dominant 
contribution (at tree-level) to the \i — > 3e decay rate and the /i — e conversion rate in a 
nucleus. In the case of the first process we have (see, e.g., [77]): 

BR(/i^3e) ~ 16|^ e | 2 ^3sm 4 e w - 2sm 2 6 w + • (93) 

Taking into account the experimental upper limit reported in (T4]), we get the following upper 
limit on the \i — e effective coupling: 

|V| < 5.6 x 10~ 7 . (94) 

which is a stronger constraint with respect to the one derived from the non-observation of 
the /i — > cy decay (see eqs. ( 1921) and (J9~TT)). mediated (at one-loop) by an effective dipole 
operator. 

More stringent constraints on the effective fi — e — Z coupling can be obtained using the 
data from the /i — e conversion experiments. Indeed, according to the general parametrisation 
given in [75] (see also [791 EH]), we have for the (/, — e conversion ratio in a nucleus Af with 
N neutrons and Z protons: 

CR(fiAf eAf) = ^ |C Me | 2 \(2g LV{u) + g LV(d) ) V® + (g LV{ u) +2g LV{d) ) V^f , (95) 
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where in this case 

C^e = 4 Vfie , (96) 

V {n) * jV^\ g LV{u) = 1-^4 and g LV[d) = + (97) 

An upper bound on \rjfj, e \ can be derived from the present experimental upper limit on the 
fi — e conversion rate in the nucleus of ^Ti, CR(/i Ti — > e Ti) < 4.3 x 10~ 12 . From eqs. ( 1951) - 
(1971) we get: 

|^ e | < 2.6 x HT 7 . (98) 

If in the fi — e conversion experiments with ||Ti the prospective sensitivity to CR(/i Ti — >■ 
eTi) ~ 10~ 18 will be reached, these experiments will be able to probe values of |?7 M e| as small 
as |^ e | ~ 1.3 x KT 10 . 

5 Discussion and Conclusions 

We have performed a detailed analysis of charged lepton flavour violating (LFV) processes 
— fi — > e7, fi — > 3e and fi — e conversion in nuclei — in the context of see-saw type extensions 
of the Standard Model, in which the scale of new physics A is taken in the TeV range, 
A ~ (100 — 1000) GeV. In this class of models, an effective Majorana mass term for the 
light left-handed active neutrinos is generated after electroweak symmetry breaking due 
to the decoupling of additional "heavy" scalar and/or fermion representations. We have 
analyzed in full generality the phenomenology of the three different and well-known (see- 
saw) mechanisms of neutrino mass generation, in their minimal formulation: i) type I see-saw 
models, in which the new particle content consist of at least 2 RH neutrinos, which are not 
charged under the SM gauge group; ii) type III see-saw models, where the RH neutrinos are 
taken in the adjoint representation of SU{2) L with zero hypercharge; Hi) type II see-saw (or 
Higgs triplet) models, where the scalar sector of the theory is extended with the addition of 
at least one scalar triplet of SU (2) L coupled to charged leptons. 

Under certain conditions the couplings of the SM charged leptons with the new fermions 
and/or scalars are, in principle, sizable enough to allow for their production and detection 
at present collider facilities, LHC included. On the other hand, remarkable indirect tests 
of such scenarios are also possible in ongoing and future experiments looking for charged 
lepton flavour violation. Indeed, the flavour structure of the interactions between the SM 
leptons and the new "heavy" particle states is mainly determined by the requirement of 
reproducing neutrino oscillation data, in such a way that the unknown parameter space can 
be expressed in terms of very few quantities. The latter can, therefore, be constrained by 
the measurement of different LFV observables. Further and complementary information is 
provided also by experiments which search for lepton number violating phenomena, such as 
neutrinoless double beta decays of even-even nuclei. 

We summarize below the phenomenological implications of a possible observation of the 
LFV processes given above for each kind of (TeV scale) see-saw extensions of the SM. 

13 The expression for V( n ) is valid under the assumption of equal proton and neutron number densities in 
the given nucleus [75] . The numerical value of the nuclear form factors for ^Ti, f 3AI and 7g 7 Au is reported 
in Tabled 
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Type I see-saw results. In this case the /x — > cy and [i — > 3e decay branching ratios 
BR(/i y ej) and BR(/i — > 3e), and the fi—e conversion rate in a nucleus Af, CR(/i A/" — > e A/"), 
A/" = Al, Ti, Au, can have values close to the existing upper limits and within the sensitivity 
of the ongoing MEG experiment searching for the /j, — > cy decay and the future planned fi — e 
conversion and /i — > 3e decay experiments [2H1 EHJ ED E21 E3] • The relevant LFV observable 
in the minimal scenario, with the addition of only two RH neutrinos to the SM particle 
content, is provided by the quantity \(RV)* fll (RV) e i\, where (RV)ij (j = 1,2) denote the 
couplings of the fermion singlets to the SM charged leptons (see eqs. f lT5|) and fTTB"]) ). If the 
MEG experiment reaches the projected sensitivity and no positive signal will be observed 
implying that BR(/i — > cy) < 10~ 13 , there still will be a relatively large interval of values 
of \(RV)*^ l (RV)ei\, as Fig. |2] shows, for which the /x — e conversion and /x — > 3e decay are 
predicted to have observable rates in the planned next generation of experiments. 

It follows from the analysis performed by us that as a consequence of an accidental 
cancellation between the contributions due to the different one loop diagrams in the /x — e 
conversion amplitude, the rate of /x — e conversion in Al and Ti or in Au can be strongly 
suppressed for certain values of the see-saw scale M 1 . As we have seen, this suppression can 
be efficient either for the conversion in Al and Ti or for the conversion in Au, but not for all 
the three nuclei, the reason being that the values of M\ for which it happens in Al and Ti 
differ significantly from those for which it occurs in Au. In both the cases of Al or Ti and 
Au, the suppression can be effective only for values of M\ lying in very narrow intervals (see 
Figs. ElandlU). 

In the case of IH light neutrino mass spectrum, all the three LFV observables, BR(/i — > 
ej), BR(/i — > 3e) and CR(/iA/" — > eAf), can be strongly suppressed due to the fact that the 
LFV factor {(RV)^ 2 oc {U^ + i^/m-i/^U^ 1 2 — WijZ + iU^ 2 , in the expressions of the 
three rates can be exceedingly small. This requires a special relation between the Dirac and 
the Majorana CPV phases 5 and «2i, as well as between the neutrino mixing angle #13 and 
the phase 5 (see eq. fl34"l) ). For the values of sin 6*13 from the current 3a allowed interval, 
eq. ([I]), one can have [U^ + iU^ 2 = provided < 5 ^ 0.7. A priori it is not clear 
why the relations between 5 and 0^21, and between 5 and #13, leading to \U^a + iU^i\ 2 = 0, 
should take place (although, in general, it might be a consequence of the existence of an 
approximate symmetry). The suppression under discussion cannot hold if, for instance, it 
is experimentally established that 5 is definitely bigger than 1.0. That would be the case if 
the existing indications [9] that cos 5 < receive unambiguous confirmation. 

We note finally that for M l ^ 100 GeV we have: BR(/i 3e)/BR(/x erf) ^ 0.067. 
Thus, if it is experimentally established that BR(/i — > 3e) /BR(/i — > e^y) is much smaller than 
the quoted lower bound, the model considered with Mi ^ 100 GeV will be ruled out. Such 
a result would be consistent also just with a see-saw scale M 1 < 100 GeV. 

Type II see-saw results. It follows from the results obtained in Section [3] that the 
predictions for the \i — > ey and /i — > 3e decay branching ratios, as well as the \i — e conversion 
rate in a nucleus A/", in the TeV scale type II see-saw scenario considered exhibit, in general, 
different dependence on the masses of the singly- and doubly-charged Higgs particles A + and 
A ++ , which mediate (to leading order) the three processes. For = m&++ — Ma, all 
the three rates are proportional to M^ 4 , i.e., they diminish as the 4th power of the see-saw 
scale when the latter increases. 

The matrix of Yukawa couplings hu> which are responsible for the LFV processes of 
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interest, is directly related to the neutrino Majorana mass matrix and thus to the PMNS 
neutrino mixing matrix U. As a consequence, BR(/i — > ej), BR(/i — > 3e) and CR(fiM — > 
eM) depend, in general, on the neutrino mass and mixing parameters, including the CPV 
phases in U . 

To be more specific, BR(/i — > ej) does not depend on the Majorana CPV phases and 
on min(mj), and its dependence on the Dirac CPV phase and on the type of neutrino mass 
spectrum is insignificant. In contrast, both BR(/i — > 3e) and CK(fiM —> eM) exhibit 
very strong dependence on the type of neutrino mass spectrum and on the values of the 
Majorana and Dirac CPV phases. As a consequence, the predictions for BR(/i — > 3e) and 
CR(/iA/" — > eM) for given Ma can vary by several orders of magnitude not only when the 
spectrum changes from NH (IH) to QD as a function of the lightest neutrino mass, but 
also when one varies only the values of the CPV phases keeping the type of the neutrino 
mass spectrum fixed. All the three observables under discussion can have values within the 
sensitivity of the currently running MEG experiment on the fi — > e'y decay and the planned 
future experiments on the n — > 3e decay and fi — e conversion. However, for a given see-saw 
scale in the range of ~ (100 — 1000) GeV, the planned experiments on /i — e conversion in Al 
or Ti will provide the most sensitive probe of the LFV Yukawa couplings of the TeV scale 
type II see-saw model. 

Type III see-saw results. Unlike the type I see-saw extension of the SM discussed in 
Section [21 in this scenario we have several — possibly sizable — lepton flavour violating in- 
teractions in the low energy effective Lagrangian, due to the higher SU(2)l representation 
of the new fermion fields. In particular, FCNCs arise at tree-level from the non-unitarity of 
the PMNS matrix (see eq. ( 182|) ). Thus, the effective fi — e — Z coupling in ( 182]) makes it 
possible an enhancement of at least two orders of magnitude of the rates of fi — > e'y, fi — > 3e 
and fi — e conversion with respect to the ones predicted in the type I see-saw scenario, with 
RH neutrinos taken in the TeV range. Consequently, all the predicted LFV observables may 
be probed in the related present and future experiments. As in the previous scenarios, the 
strongest constraint on the flavour structure of this class of models is by far provided by the 
expected very high sensitivity reach of fi — e conversion experiments. 

In conclusion, the oncoming combination of data on neutrino oscillations, collider searches 
and lepton number /flavour violating processes represent an important opportunity to reveal 
in the next future the fundamental mechanism at the basis of the generation of neutrino 
masses as well as the underlying physics beyond the standard theory. 
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